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With current magnetic field shielding and high precision detection in dipolar spinor Bose-Einstein 
condensates, it is possible to experimentally detect the low or zero field nonsecular dipolar dynam¬ 
ics. Here we analytically investigate the zero-field nonsecular magnetic dipolar interaction effect, 
with an emphasis on magnetization dynamics in a spin-1 Bose-Einstein condensate under the sin¬ 
gle spatial mode approximation within the mean field theory. Due to the biaxial nature of the 
dipolar interaction, a novel resonance occurs in the condensate magnetization oscillation, contrast 
to the previous assumption of a conserved magnetization in strong magnetic fields. Furthermore, 
we propose a dynamical-decoupling detection method for such a resonance, which cancels the stray 
magnetic fields in experiments but restores the magnetization dynamics. Our results shed new lights 
on the dipolar systems and may find potential applications beyond cold atoms. 

PACS numbers: 03.75.Kk, 03.75.Mn, 67.85.Fg 


I. INTRODUCTION 

Magnetic dipolar interactions exist in a wide variety of 
physical systems, including nuclear spins [l|, atoms and 
molecules ii , condensed matter Q, chemical @ and 
biological systems 0. Due to the relative weakness of 
the magnetic dipolar interaction to the Zeeman effect of 
an external magnetic field, the nonsecular interaction in 
a dipolar system, which varies temporally with Larmor 
frequency and breaks the rotational invariance, was usu¬ 
ally neglected in most experiments over decades. Such 
an approximation results in the famous magnetization 
conservation during the time evolution [I|. 

Recently, with the precision measurement of an ultra- 
weak or zero magnetic field, the investigation of the non¬ 
secular dipolar interaction effects becomes revived in the 
areas of zero-field nuclear magnetic resonance and dipo¬ 
lar spinor Bose-Einstein condensates [7l-[TTj| . Particularly, 
the achievements in spinor Bose-Einstein condensates 
(BECs) provide a highly tunable and controllable sys¬ 
tem where the spin interactions, including the magnetic 

S olar interaction, can be accurately engineered SIH- 
. Such dipolar spinor condensates make them an ideal 
testbed to extensively explore the full magnetic dipolar 
interaction effect, including not only the secular part but 
also the nonsecular one which breaks the rotational sym¬ 
metry and definitely brings new physics. Recent experi¬ 
ments by Pasquiou et al. investigated the incoherent de¬ 
magnetization process i n sp in-3 ®^Cr condensates in an 
ultralow magnetic field [1^. However, the coherent dy¬ 
namics induced by the nonsecular dipolar terms in cold 
atoms has yet been touched. 

In this paper, we investigate the effect of the magnetic 
dipolar interaction in a spin-1 BEG, focusing on the co¬ 
herent magnetization dynamics induced by the nonsecu¬ 
lar dipolar interaction in low or zero field (See Fig. [2]). 


By adopting mean field theory and single spatial mode 
approximation (SMA) [20l - l^ . we analytically obtain the 
magnetization dynamics of the spin-1 BEG. Gontrast to 
the usual assumption of magnetization conservation, we 
find a novel resonance in the magnetization oscillation 
in the dipolar condensate. With large-scale numerical 
calculations, we also explore the spin dynamics beyond 
the SMA, which agrees well with the analytical predic¬ 
tion. Furthermore, a practical control protocol utilizing 
dynamical decoupling techniques is proposed to liberate 
the nonsecular part of the magnetic dipolar interaction 
by canceling the stray magnetic fields. Our results shed 
new light on the dipolar systems, in particular on the 
resonant magnetization dynamics induced by the non¬ 
secular part of the dipolar interaction, and open a door 
to understanding the complex spin texture observed in 
dipolar spinor BEGs |23l - l2^ . 

The paper is organized as follows. In Sec. |TT1 we intro¬ 
duce the mean field Hamiltonian describing the dipolar 
spin-1 condensate at zero magnetic field. In Sec. IHIl 
we present the equation of motion for the condensate 
magnetization and the analytical solution. We then de¬ 
sign dynamical decoupling pulse sequence to suppress the 
Zeeman effect of a nonzero magnetic field and perform 
numerical simulation with parameters in practical exper¬ 
iments beyond SMA in Secs. IIVI and IVl respectively. We 
conclude in Sec. lYD 


II. SYSTEM HAMILTONIAN 

In the mean field theory and under the SMA, a dipo¬ 
lar spin-1 Bose condensate is described by the Hamilto¬ 
nian 

H = cf + ds{3f! - f) + 3dMy - /J), 
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where the spin exchange interaction strength is c = 
{C 2 I 2 ) f drp^(r) with p(r) = |</>(r)p being the density 
of the mode function ^(r). The spin-exchange interac¬ 
tion coefficient C 2 = 47r/i^(a2 — ao)/(3M) with M being 
the atom’s mass and 09,2 the s-wave scattering length 
of two spin-1 atoms in the symmetric channel of the 
total spin 0 and 2, respectively. The condensate spin 
P = fx+fy+ /z^ where = P\F^,y,z\^ with |f) = 

being the spin wave function and F^^y^z the 
spin-1 matrices. Note that we have neglected the con¬ 
stant terms under the SMA, contributed by the kinetic 
energy, the trapping potential, and the spin-independent 
interaction proportional to cq = 47r/i^(2a2 -I- ao)/(3M). 

The dipolar interaction includes two parts, the secu¬ 
lar part with dg = (cd/4) / drdr'jr — r'|“^p(r)p(r')(l — 
3 cos^ Oe) and the nonsecular part with dn = 
(C(i/4) / drdr'jr — r'|“^p(r)p(r') sin^ where 6e 

and ipe are respectively the polar and azimuthal an¬ 
gles of (r — r'). The dipolar interaction coefficient is 
Cd = MoMB 5 F/( 47 r) with po being the vacuum permeabil¬ 
ity, pb the Bohr magneton, and gp the Lande g-factor 
for an electron. For ®'^Rb atoms, Cd/|c 2 | ~ 0.09, and 
for ^^Na atoms, Cdl\c 2 \ ~ 0.007. Note that = 0 if 
p(r) is axially symmetric. In general, the value of dg.n 
can be positive, 0 , or negative, depending on the specific 
shape of the mode function [^, . In a large or mod¬ 

erate magnetic field, the secular part survives due to the 
conservation of /z, but the nonsecular part, which flips 
the spin, is strongly suppressed by the Zeeman effect and 
usually neglected. In Sec. IIVI we will propose a dynam¬ 
ical decoupling method to cancel the Zeeman effect and 
to “restore” the suppressed nonsecular part. 

The original Hamiltonian can be rewritten as 



FIG. 1. (Color online) The iso-energy contour plot of the 
anisotropic energy in unit of |c| on the surface of a spin 
sphere (a) for category (I) with Aa, = 0.5 and Az = 1 and (b) 
for category (II) with Aj, = — 1 and A 2 = 1. The black lines 
mark the boundary between the « and y regions in (a) and 
between the 2 and x regions in (b). 

nothing but a nonlinear rigid pendulum, which oscillates 
around x-, y-, or z- axis correspondingly. 

III. ANALYTICAL RESULTS OF 
MAGNETIZATION RESONANCE UNDER SMA 

By treating the condensate spin as a classical spin, 
which rotates in an effective magnetic field {bx,by,bz) = 
{2Axfx,0,‘^Azfz), we obtain the following equation of 
motion for the magnetization, = df^/dt = 2Axfxfy 
By utilizing again the relations fy = p—p—fz P = 
{Fxz—Azp)IAj;, we find a closed equation of motion for 
fy 

P = A{Azp - Exz) [(A, - Az)p + Exz - Axp] , (2) 


H = AoP+Axp + Azp (1) 

where Aq = c—ds+3dn, A^ = —6dn, and A^ = 3{ds—dn)- 
We have used the relation that p = p—p—p- For such 
a spin system, besides the conservation of the total energy 
and the spin / of the condensate (thus the conservation 
of the isotropic energy Eq = AqP), we notice that the 
sum of the rest two anisotropic terms, which we define as 
Exz = ^xP + ^zP, is also conserved. The conservation 
of Exz defines the spin trajectories in the x-z plane into 
the following two categories, depending on the values of 
Ax and A^ (besides the trivial one, AxAz = 0 where fz 
or fx is constant): (I) Aa;A 2 > 0, the trajectory in the 
x-z plane is an ellipse (or a circle for Ax = Az); (II) 
AxAz < 0, the trajectory is a hyperbola (or a parabola 
for Ax = —Az). Of course, fx,z niust be smaller than 
{fy 0 ) or equal to / {fy = 0 ). 

The iso-energy contour plots of Exz on the surface of 
the spin sphere for categories (I) and (II) are presented 
in Fig. [T] For both categories the oscillations of fz either 
passes zero or not, depending on the initial value of /^. 
As shown in the figure, there is an alternative way to 
classify the spin dynamics: the x-, y-, or z-region, where 
the spins rotates around the corresponding x-, y-, or z- 
axis. Actually, the trajectory of the spin behaves like 


where / and Exz are determined by the initial state. 

The analytical solution for fz is an inverse function of 
the elliptic integral of the first kind A[-, •], 


t — to + 


2pAz{Axp-Exz) 


Sin 




(3) 


where x'^ = Exz!Az and = {AxP - Exz) j{Ax - A^). 
The value of x (y) is determined by fx{y) = 0 (i.e., fz = 
0 ). 

The period of the oscillation is 


T= ifz{t)-^df, 
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FIG. 2. (Color online) Periodic spin dynamics induced by the 
magnetic dipolar interaction. Oscillation of the magnetization 
fz with an initial rotation angle (a) 0 — 15° (below the critical 
value 6c) and (d) 9 = 75° (above 6c), where the spin oscillates 
around a- and x-axis, respectively, (b) Dependence of the 
oscillation period T on the initial angle 6. Solid lines are 
the SMA prediction and circles are the numerical results with 
coupled Gross-Pitaevskii equations. A resonance occurs at 6c, 
marked by the black arrow in (b), where the spin evolves along 
neither x- nor 2 - axis, (c) Typical iso-energy spin trajectories 
of the condensate on the surface of the spin sphere. The red 
and blue arrows show the spin evolution direction. 


if in the z-region, where K[-) is the complete elliptic in¬ 
tegral of the first kind. The function K{k) is nearly a 
constant 7r/2 around k = 0 and diverges rapidly if fc —>■ 1. 
In addition, K{k) = F{'k/ 2, k), and K{\lk)lk^l'^ also di¬ 
verges if fc —>■ 0. 

The analytical solution for shows oscillatory and 
periodic motion. This is in sharp contrast to the previous 
understanding of the conservation of the magnetization 
(in large magnetic fields), because the nonsecular terms 
break the rotational symmetry around z-axis. As shown 
in Fig.[Tl the amplitude of the oscillation of is |?/| if the 
condensate spin is in the a:-region, |a:| if in the y-region, 
or \x — y\ if in the z-region. Remarkably, the period T 
diverges if a; = 0 (j/ = 0), i.e., the initial spin state is 
set on the boundary between the x-region (y-region) and 
the z-region (see Fig.[T]). This divergence indicates that a 
resonance occurs by changing the initial spin state across 
the boundary. Such an interesting resonance has not been 
revealed before, because of the neglect of the nonsecular 
terms, and is obviously due to the competition between 
the secular and nonsecular dipolar terms. 

For a clear view, we present in Fig. [2] the analytical 
results for the category (II) with ds/\c\ = 0.225 and 
dn/\c\ = 0.0644 (|c| = 1.3 Hz is calculated for a numeri¬ 
cally simulated spin-1 BEG). Two typical oscillations of 
fz in the z-region and x-region are shown in Fig. HJa) 
and (d). Other trajectories are shown in Fig. [2Kc) and 
the periods of the oscillations with respect to initial polar 


angle 6 {f = 1 and zero azimuthal angle) of the conden¬ 
sate spin are shown in Fig. [2Ib). We observe a clear res¬ 
onance signature at 9c, which is determined by E^z = 0, 
i.e., cos^dc = Axl{Ax — Az). Similarly, for the category 
(I) we also observe a resonance in the fz oscillation period 
and the critical angle cos^ 9c = A^/Az, which is obtained 
from Az:f^ - E,zz = 0 [3l|. 

The resonance in the oscillation period can be under¬ 
stood physically by treating the dipolar spin-1 BEG as 
a two-axis rigid nonlinear pendulum (see Eig. [5]). For 
a single-axis pendulum, i.e., either A^ = 0 or = 0, 
the condensate spin oscillates around an axis, z or x, 
with a constant spin component along this axis. But 
in categories (I) and (B), neither of A^^z is zero, so the 
evolution of the condensate spin is not purely around a 
single axis but a more complicated oscillation between 
the two axes, due to the competition of the two nonlin¬ 
ear terms of A^f^ and Azff. In general, by setting the 
condensate spin initially closer to an axis, e.g., z-axis, the 
evolution of the condensate spin is around this axis since 
the corresponding term (Azff) is dominant. However, 
there exists a clear boundary where the two-axis terms 
are balanced, e.g, A^ff = —Azff in Fig. HKb), and the 
condensate spin evolves in a third direction and reaches 
to a dead end {fx,z = 0). The oscillation never completes 
and the period becomes infinite. Thus a resonance peak 
appears in the spin oscillation period. 

An alternative way of understanding the resonance 
behavior is to treat the condensate in a semiclassical 
manner. In a single particle picture, the energy level 
of an atom in the {mp = ±1) state is shifted by an 
amount of Az due to the term Azff. Similarly, the 
term A^f'^ not only shifts oppositely the energy level 
of |mF = ±1) (if AxAz < 0), but also couples them. As 
shown in Fig. [2] and [TJb), when the resonant condition 
Axfx + ^z/z = 0 is satisfied and the condensate reaches 
its steady state {fx,z = 0), the total shift of the energy 
levels of \mp = ±1) become zero, i.e., the energy levels 
of ImF = ±1) are resonant to that of ImF = 0). 


IV. SUPPRESSING THE ZEEMAN EFFECT OF 

AN EXTERNAL STRONG MAGNETIC FIELD 

Inside a magnetic shield room, the stray magnetic field 
have been reduced to as low as O.I inG [ifljllll- However, 
even in such a low magnetic field, the Zeeman effect still 
overwhelms the dipolar-interaction-induced spin dynam¬ 
ics along the z-axis, since the Zeeman energy /isB ^ 140 
Hz for B = 0.1 mG is much larger than the dipolar inter¬ 
action energy, typically in the order of 0.1 Hz at a con¬ 
densate density of 10^"^ cm“^. As shown in Figs. EJa), 
the oscillation amplitude of /^ for B = 0.1 mG (dash- 
dotted line and the inset) is much smaller than that for 
B = 0 (dotted line). The consequence is that the magne¬ 
tization dynamics due to the dipolar interaction is rather 
challenging to observe in a practical experiment because 
the stray magnetic held effects dominate the magnetiza- 
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FIG. 3. (a) Coherent magnetization dynamics of a dipolar 
spin-1 BEG in a magnetic field. B = 0 (dotted line), B = 0.1 
mG (dash-dotted line), B = 0.1 mG with control pulse delay 
r = 0.004 s (solid line), and B = 0.1 mG with r = 0.04 s 
(dashed line). The initial condition is the same as in Fig. [2] 
with 9 = 30°. Other parameters are the same as in Fig. (2] 
The Zeeman effect of the nonzero magnetic field suppresses 
the fz dynamics (dash-dotted line and the inset for a zoom- 
in view), but the fast and periodic application of ttx. pulses 
restores the dipolar fz dynamics, (b) Dependence of the nor¬ 
malized amplitude of the dipolar magnetization oscillations 
on the control pulse delay r at B = 0.1 mG (solid line) and 
B — 1 mG (dashed line). 


tion dynamics. 

However, the resonant dipolar magnetization dynam¬ 
ics can be revealed under current experimental conditions 
if we employ dynamical decoupling techniques to cancel 
the Zeeman effect [U, In particular, by applying 

frequent TTa, pulse, which rotates the condensate spin 180° 
along a;-axis, the Zeeman effect is eliminated while leav¬ 
ing the dipolar interaction intact to the leading order. 
We assume that the control pulse is an instantaneous 
(hard) pulse and the delay between two adjacent pulses 
is T (Appendix O . The numerical simulation results are 
presented in Fig. |31 As shown in Fig. [SJa), in the limit 
of small T (solid line), the dipolar magnetization dynam¬ 
ics is restored; in the limit of large r (dashed line), the 
Zeeman effect is not well suppressed and the dipolar mag¬ 
netization oscillation amplitude is small. Under control 
pulses, the dependence of the normalized magnetization 
oscillation (which is the ratio of the magnetization oscil¬ 
lation amplitude at H = 0.1 mG under control pulses to 
the free magnetization oscillation at H = 0) on the pulse 


delay r is shown in Fig. [3Kb). Obviously, the Zeeman ef¬ 
fect is canceled and the magnetization oscillation is well 
restored at i? = 0.1 mG if r is smaller than 0.02 seconds, 
which is easily realizable in experiments (Appendix lAl). 


V. NUMERICAL RESULTS WITH A LARGE 
NUMBER OE ATOMS 

Our previous analysis are based on the SMA, which 
may not be valid in some experiments for spin-1 BEGs, 
particularly for a condensate with a large number of 
atoms [13, [H, 113 . It is an open question whether 
the SMA results remain valid under practical experi¬ 
mental conditions, where the number of ®^Rb atoms 
is 10^ and the trap frequencies are {wa,,Wy,Wz} = 
(27rx){90,140, 200} Hz. For this trap geometry, the SMA 
with a three-dimension Gaussian wave function predicts 
that ds = 0.215 and dn = 0.066, while the numerical re¬ 
sults for a fully polarized ground state (all atoms are in 
I -b 1) state) are dg = 0.225 and = 0.064. It shows 
obviously that the SMA is slightly violated. 

We start the numerical simulation by tilting the fully 
polarized ground state away from the z-axis by a polar 
angle 0 in the x-z plane. In this way, the subsequent 
magnetization evolution is solely due to the dipolar in¬ 
teraction since the condensate magnetization would re¬ 
main constant in the absence of the anisotropic dipolar 
interaction. The magnetization dynamics of the conden¬ 
sate is obtained by numerically solving the three coupled 
Gross-Pitaevskii equations in the given trap [13,113, HE] 

where the kinetic energy is T = —/i^V^/(2M) with 
M being the atom mass of ®^Rb, the trapping poten¬ 
tial is Uext(a;,?/, z) = M{ujIx'^ + and 

the total number density is n = Ha V'aV’a with 
(a = ±1,0) being the three components of the con¬ 
densate wave function. The effective field originat¬ 
ing from the spin-exchange and dipolar interactions is 
Beff = C2S-|-Cd/dr'{S(r')- 3 [S(r')-e]e}/|r-r'p, where 
S = ip^Fa/sipiS is the spin density with F the atom spin-1 
matrix, and e is the unit vector along r — r'. The spin- 
exchange interaction coefficient C2 and the dipolar inter¬ 
action coefficient Cd are given previously. We solve these 
coupled equations using the operator splitting method, 
where the term involving the integral operator Beff is cal¬ 
culated with convolution theorem and fast Fourier trans¬ 
form. 

We present in Figs. UKa), HKb), and HKd) the numeri¬ 
cal results, as well as the SMA predictions with dg^n ob¬ 
tained numerically from the fully polarized ground state. 
We observe a pretty good agreement of the spin oscilla¬ 
tions between the numerical results and the SMA predic¬ 
tion, due to the fact that the three components of the 
ground state of the ferromagnetically interactin g sp in-1 
condensate share the same spatial wave function [2lL l22j| . 
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However, we find slight mismatches of the oscillation am¬ 
plitude in Fig. [D^a) and the resonance peak position in 
Fig. EJb). These discrepancies might be due to the trap 
anisotropy: the atoms spins are more likely to align along 
the loosely trapped x- or y-axis so that the total dipolar 
energy is lower, which means that the z-region is smaller 
than the SMA results. Consequently, the boundary be¬ 
tween the a:-region and the z-region is shifted to a smaller 
6c, which is what we observe in Fig. [^b), and the lower 
value of fz in Fig. ^a) becomes smaller. 


VI. CONCLUSION 

The weak dipolar interaction effect, particularly the 
effect of the nonsecular part, in spin-1 Bose conden¬ 
sates such as ®^Rb is experimentally challenging to ob¬ 
serve [13, . Our proposal provides a practical way to 

detect the condensate magnetization oscillation induced 
by the nonsecular dipolar interaction, where a resonance 
emerges in the oscillation period. Numerical results with 
experimental parameters, such as the large number of 
atoms and the suppression of Zeeman effect of the stray 
magnetic field with dynamical decoupling method, con¬ 
firm that the resonant behavior of the coherent magne¬ 
tization dynamics is detectable experimentally. Our re¬ 
sults point to a new direction for future investigations 
on many dipolar effects, including the competition with 
the short-range contact interaction, the quantum dipolar 
effects which are of great interests for quantum metrol¬ 
ogy and next-generation-precision magnetometers based 
on spinor Bose condensates [ll|,I13“l4l| , and the structure 
determination in chemical and biological dipolar systems. 
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Appendix A: Dynamical decoupling 

For a spin-1 condensate in a magnetic field, H = 
<-^zSz + V, where uJz is the Zeeman splitting and V is the 
weak two-boday interaction (including the spin exchange 
and dipolar interaction). Suppose the control pulse 
is an instantaneous (hard) pulse and the delay between 
two adjacent pulses is r. Note that the V is unaffected 
but Sz becomes —Sz by the tTx pulses. Let us consider 
the evolution after 2 pulses [I|, 

U{t = 2T)=U{T)U^^U{T)U^^. 

Since = exp(—i7r5x) = and U{t) = 

exp[—iT{uJzSz + V)] ~ ex^{—iTU}zSz) exp(—irU) for 
small r (i.e., ||iL||T ^ 1), we find that 

U{t = 2T)=U{T)U^^U{T)Ul 

Ki U(T)exp(-iTU)t/^,„ ex-p{-iTWzSz)Ul^ 

~ exp(—i2TU). (Al) 

After periodically applying 2K pulses, the evolution 
operator is 

U{t = 2 Kt) = [C/(T)[7^j2ic 

~ exp{—itV) (A2) 

to the leading order of r. Clearly, the effect of the mag¬ 
netic field is removed if t —> 0. 

A possible way to realize experimentally the hard pulse 
is by applying a square pulse of constant Bi, e.g. 10 mG, 
along X direction within a pulse width w ~ 10“^ s. The 
pulse width is much smaller than the typical pulse delay 
T 10“^ s and the dipolar dynamics time scale ^ 10 

sna. 
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